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OSTEOPOROSIS DIAGNOSIS AIDING APPRATUS UTILIZING PANORAMIC 

RADIOGRAPHS 

Technical Field 

5 The present invention relates to an apparatus for carrying out diagnostic aid of 

osteoporosis utilizing panoramic radiographs taken during dental care. 

Background Art 

Early detection and early treatment of patients suffering from osteoporosis 
10 estimated to be a population often million and a potential group thereof currently in Japan 
is important. 

Current major osteoporosis diagnostic techniques in Japan are- bone-density 
measurement of the lumbar spine and the femur by a double X-ray resorption 
measurement apparatus (DXA) capable of measuring bone density with utmost accuracy; 

15 and bone density measurement of the calcaneus through quantitative ultrasonography 
(QUS) as a simplified method. However, the current state is that only certain facilities 
are equipped with these apparatus, and are only used for persons with symptoms of 
osteoporosis such as lower back pain, or persons with great concern for osteoporosis, while 
not being used for prevention of fractures considered to be caused by osteoporosis. 

20 Recent studies claim that a decrease in bone density may bring about alterations 

in shape of the cortical bone at the lower edge of a mandible (thick compact bone covering 
the outer side of the mandible bone) in a panoramic radiograph of all teeth and 
surrounding bones of the chin (see Non Patent Reference 1). FIG. 1 shows an exemplary 
panoramic radiograph of teeth. In addition, it is understood that persons with an altered 

25 shape of the lower edge of the mandible have a significantly higher risk of fracture due to 
slight external pressure compared to persons without (see Non Patent Reference 2). 

Common dental care does not put diagnosis of osteoporosis using panoramic 
radiographs into practical use since examination of this mandible inferior border area is 
infrequent and an examiner's reading ability of X-rays is questioned. 

30 Although studies on diagnosing osteoporosis utilizing panoramic radiographs 

have mainly been centered on altered shapes of a portion called spongy trabecular bone 
existing in the maxillary bone and the mandible bone (see Non Patent References 3, 4, 5, 
and 6), diagnosis by automatically extracting only the trabecular bone from a panoramic 
radiograph is difficult, and therefore has not been put into practical use. It has also been 

35 reported that diagnosis is possible from the thickness of the cortical bone at the lower edge 
of the mandible (see Non Patent References 3 and 7). 



If persons at risk of osteoporosis can be screened using panoramic radiographs, 
screening of persons without subjective symptoms that are difficult to diagnose or persons 
without concern for osteoporosis, and instruction for only persons potentially having 
osteoporosis to undergo closer examination such as DXA or referral to a facility equipped 
5 with that apparatus allow early detection and early treatment of patients suffering from 
osteoporosis, and also reduce costs of the examination. Furthermore, the method used for 
this screening must be simple and usable even without having any special skills or 
requiring complicated operations. 

However, technology for diagnosing a patient suffering from osteoporosis by, for 
10 example, focusing on an altered shape of the cortical bone at the lower edge of the 
mandible in a panoramic radiograph, and objectively and automatically assessing it 
through mathematical morphologic image processing using computer software does not 
exist. 
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Disclosure of the Invention 
An objective of the present invention is to provide an osteoporosis diagnosis aiding 
apparatus utilizing panoramic radiographs which can screen persons at a risk of 
osteoporosis through simple operations. 
5 To reach the above-given objective, the present invention is an osteoporosis 

diagnosis aiding apparatus utilizing a panoramic radiograph of teeth is characterized by 
comprising: image region extracting means for extracting an image of a region including a 
specified portion in the panoramic radiograph of teeth; skeleton processing means for 
finding a skeleton of the extracted image; line extracting means for extracting only 
10 portions parallel to the tilt of the lower edge of the mandible from the found skeleton^ 
binarizing means for binarizing an image including the skeleton so that lines of the 
skeleton are separated from the background; and determining means for determining 
whether there is a risk of osteoporosis by classifying the binarized lines according to size. 

The determining means can classify the binarized lines into three groups 
15 according to size, and determine there is a risk of osteoporosis in a case where there is 
more than just a single line classified in the largest line group aside from the smallest line 
group. 

The image extracting means may extract an image of a rectangular region parallel 
to the tilt of the lower edge of the mandible, and the line extracting means may extract 
20 horizontal portions of the lines. 

A program designed to establish processing of the osteoporosis diagnosis aiding 
apparatus for a computer system also falls within the scope of the present invention. 



Brief Description of Drawings 
25 FIG. 1 is a diagram showing an exemplary panoramic radiograph of teeth; 

FIG. 2 is a diagram showing a mandible according to a panoramic radiograph, a 
normal cortical bone and a reduced bone density patient's cortical bone; 

FIG. 3 is a flowchart showing image processing of an osteoporosis diagnostic 
aiding apparatus; 

30 FIG. 4 is a diagram showing image processing for a normal person; 

FIG. 5 is a diagram showing image processing for a person with osteoporosis; 

FIG. 6 is a diagram describing the meaning of skeleton; 

FIG. 7 is a diagram describing processing for providing a skeleton; 

FIG. 8 is a graph showing a histogram of linear elements; and 
35 FIG. 9 is a graph showing a partially enlarged histogram of the linear elements. 
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Best Mode for Carrying Out the Invention 
An embodiment of the present invention is described forthwith while referencing 
the appended drawings. 

The present invention is for screening patients suffering from osteoporosis by 
5 identifying altered shapes of cortical bone at the lower edge of the mandible in panoramic 
radiographs using a skeleton, which is a mathematical morphologic image processing 
method, based on reports of Klemetti et al. (see Non Patent Reference l). 

A method according to the present invention allows screening of a person without 
subjective symptoms or concern for osteoporosis since panoramic radiographs of a patient 
10 who has visited a dental office for dental care are used, and it does not require special 
operations or skills since diagnosis thereof is carried out semrautomatically. Note that 
use of a skeleton in panoramic radiographs is described in Non Patent Reference 8. 

FIG. 2A is a schematic diagram of the mandible in the panoramic radiograph of 
FIG. 1. The inner surface of cortical bone at the lower edge of the mandibular molar of 
15 the person who does not have decreased bone density indicated by a dashed line in the 
panoramic radiograph of FIG. 2A is smooth and has a fixed thickness (see FIG. 2B). 
However, it is evident that the inner surface linearly resorbs when bone density decreases, 
making a rough structure (see FIG. 2C). Note that as shown in FIG. 2C, trabecular bone 
in the upper portion of the cortical bone also decreases in bone density, thereby altering in 
20 shape. 

The present invention focuses on such altered shape, and detects a person at risk 
of osteoporosis with decreased bone density by semrautomatically determining whether 
the inner surface of cortical bone at the lower edge of the mandibular molar is smooth or 
rough. 

25 This method allows easier diagnosis of osteoporosis as compared to the diagnosis 

methods based on the trabecular bone or the thickness of the mandibular cortical bone. 

The present invention may be implemented using a personal computer (PC) 
program. 

The method of the present invention is described based on a flowchart for 
30 processing according to the embodiment of the present invention shown in FIG. 3, and 
extracted images shown in FIGS. 4 and 5. FIG. 4 is an image of a normal person, and FIG. 
5 is an image of a person with osteoporosis. 

First, a panoramic radiograph (see FIG. l) to be used needs to be converted to a 
digital image readable by a personal computer using a scanner or the like. Recommended 
35 resolution at this time is 300 dpi or greater. Recently, the conversion using a scanner or 
the like is not necessary as long as digital imaging devices being popularized are used. In 
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this manner, a digitalized image of the panoramic radiograph enters the personal 
computer (S210). 

The cortical bone at the lower edge of the mandibular molar in the panoramic 
radiograph is specified with a mouse to make that part a target of examination (S220). 
The part to be specified is preferably a portion below the first molar tooth. As a result, an 
area of 100 x 100 pixels centered on the specified part is extracted, for example (see FIGS. 
4A and 5A). Note that the extracted images shown in FIGS. 4A and 5A result from 
finding a tangent of the lower edge of the mandible in proximity to a specified place and 
then extracting a rectangle parallel to the tilt of that tangent when extracting that area. 
This is for facilitating subsequent extraction of only components parallel to the tilt of the 
lower edge of the mandible. 

These extracted images are respectively subjected to the following image 
processing. 

(1) Subject the image to median filtering so as to reduce noises across that image. 

(2) Find a skeleton constituted by micro-structural elements (S230). 

(3) Extract only components parallel to the tilt of the lower edge of the mandible (S240). 

(4) Binarize the image using Otsu's linear discrimination method, for example (S250). 
Each processing is described in detail below. 

<Median Filtering> 

A median filter is used to reduce noises in an image. The median filter is capable 
of removing noises that considerably differ from values in a surrounding region and are 
small and scatter randomly. This filter is known as a function of a typical image 
processing program and is used according to needs. 
<Skeleton Processing (S230)> 

Skeleton means 'bone structure*, and is to trim away the object in the image, 
creating a bone structure. Assuming the object as X, a skeleton SK (X, B) according to a 
structural element B is defined as- 
(Equation l) 
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means 'a set of centers of similar figures nB to the structural element when deploying nBs 
throughout a figure X\ 

The meaning of a skeleton is described with reference to FIG. 6 (figure X* 
rectangle, B: circle). 
5 Assuming n'<n, 
(Term 3) 

(X QnB) — {X QnB) B 

The structural element nB (left side) deployed in the center has an uncovered portion even 
if the larger n'B is deployed inside X. nB (right side) deployed in any other place can be 
10 completely covered by n'B that is inside X. 
In this manner, 

(Term 4) 

{XenB)-(XQnB) B 

has a property of 'not being able to completely cover when nB is deployed in this position 
15 even if a similar figure with a larger n is deployed inside X\ Thus, considering 
'completely covering X (e.g., a rectangle) using a similar largest figure to B (e.g., a circle)', 
nB deployed at 
(Term 5) 

{XenB)-iXQnB) B 

20 is absolutely imperative and cannot be replaced with a larger similar figure. A collection 
thereof becomes a skeleton. 

Therefore, the skeleton of figure X (rectangle) in FIG. 7A is 'a locus of the centers 
of similar largest figures to the B (circle) when completely covering X (rectangle) (see FIG. 
7B) with those similar largest figures' (see FIG. 7C). The extracted image is constituted 
25 by one or more linear elements through this image processing. 

FIGS. 4B and 5B are results of subjecting the extracted images of FIGS. 4A and 
5A to skeleton processing. 

<Extraction of Lines Parallel to Lower Edge of Mandible (S240)> 

As a result of carrying out skeleton processing, as shown in FIGS. 4B and 5B, 
30 images with fine lines running in various directions can be obtained. 

There are bundles of blood vessels in the cortical bone at the lower edge of the 
mandible. A person with decreased bone marrow may have linear bone resorption along 
the length of the luminal tissue thereof (see FIG. 2C). Since this resorption appears as a 
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line running parallel to the lower edge of the mandible, extraction of lines almost parallel 

to the lower edge of the mandible is carried out to detect this line. 

More specifically, a rectangular (linear) pattern running parallel to the lower edge 

of the mandible is prepared, and lines corresponding to that pattern are extracted. 
5 To extract the skeleton images of FIGS. 4B and 5B as rectangles running parallel 

to the lower edge of the mandible, a horizontal pattern is prepared and lines corresponding 

to that pattern are extracted. In the case of a skeleton, preparing a diagonal pattern 

running parallel to the lower edge of the mandible is difficult since the pattern for 

extraction is small and detailed. Instead, a horizontal pattern should be prepared, and 
10 image extraction is carried out in advance so as to provide rectangles running parallel to 

the lower edge of the mandible. 

<Image Binarization (S250)> 

While image processing has been carried out thus far using gray-scale images, 

now, binarization is carried out so that the lines developed through skeleton processing can 
15 be separated from the background. Binarization should be carried out by specifying a 

binarization threshold using Otsu's binarization method (see N. Otsu, 'Automatic 

threshold selection based on discrimination and least square method', IEICE Transaction, 

Vol. J63-D No.4 1980, April), for example. 

Resulting mages from extracting only portions parallel to the lower edge of the 
20 mandible from the skeleton images of FIGS. 4B and 5B and then subjecting them to 

binarization are shown in FIGS. 4C and 5C. 

discrimination of Osteoporosis (S260)> 

Linear elements appearing after binarization may be divided into three categories 

according to the number of pixels (size) thereof. 
25 In order from small to large- 

1. Noise Group 2 

2. Somewhat small linear elements that have appeared due to linear bone resorption 

Group 2 

3. Principal line of cortical bone at the lower edge of the mandible Group 3 

30 In these classifications, a threshold in number of pixels is necessary for each category. 

A panoramic radiograph is taken for a patient who has come to the department of 
radiology at the Hiroshima University Hospital dental clinic for bone-density 
measurement through DXA, and using 100 copies of the resulting image, the size of the 
linear elements after image processing is measured, creating a histogram as shown in FIG. 

35 8. As can be understood from FIG. 8, linear elements of a small number of pixels 
(approximately 0 to 100 pixels) considered to be noises are remarkably increasing. Linear 



7 



elements of 117 or less pixels or values at the base of the unimodal distribution in the 
histogram are categorized as noises (Group l). 

Next, portions less than 117 pixels are deleted, and a histogram having just larger 
components is created. The created histogram is shown in FIG. 9. FIG. 9 shows a 
5 bimodal distribution and shows that the value in between two peaks is 231 pixels. 
Accordingly, criteria that linear elements having 231 pixels or greater are categorized as 
largest linear elements of the cortical bone at the lower edge of the mandible (Group 3) and 
that linear elements having 118 to 230 pixels are categorized as those that have appeared 
due to linear bone resorption (Group 2) are made. 
10 The following judgments are made from the results of these classifications (S260). 

• The noises of Group 1 are not of concern and are disregarded. 

• Existence of linear elements of Group 2 means a risk of osteoporosis. 

• Existence of no linear elements of Group 2 and only one linear element of Group 3 means 
normal. 

15 • Existence of two or more linear elements of Group 3 means a risk of osteoporosis since a 
line due to bone resorption (line of Group 2) probably extends to be a linear element from 
Group 3. 

If the above-given judgments say that there is only a single linear element 
(principal line) categorized in Group 3 and thus diagnose as normal, 'normal' is output 
20 (S270). Otherwise, if they diagnose as a risk of osteoporosis, 'risk' is output and displayed 
on a display unit or the like (S280). 

With the examples of FIGS. 4C and 5C, FIG. 4C, since there is only a single 
principal line in the images aside from the noises, judgment of normal is output. However, 
since FIG. 5C has many long lines, judgment of a risk of osteoporosis is output. 
25 <Examples> 

Panoramic radiographs taken for 100 female patients over 50 years of age who 
came to the department of radiology at the Hiroshima University Hospital dental clinic for 
bone-density measurement of the lumbar spine and the femur through DXA are diagnosed 
using the apparatus according to the present invention. Accuracy of diagnosis by the 

30 apparatus according to the present invention is measured using these diagnostic results 
and diagnostic results through DXA according to diagnostic criteria of the Japanese 
Society for Bone and Mineral Research (persons with young average bone density under 
80% are defined as persons with decreased bone density) as gold standards. 

The results are- 82.2% sensitivity (rate of persons with a risk of osteoporosis 

35 detected by the apparatus according to the present invention to all the persons with 
decreased bone density diagnosed using DXA); 60.0% specificity (rate of persons detected 
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by the apparatus according to the present invention as having normal bone density to all 
the persons with normal bone density diagnosed using DXA); and 70.0% correct diagnosis 
rate (matching rate between all diagnostic results using DXA and all diagnostic results by 
the apparatus according to the present invention.) 
5 Twenty-seven practicing dentists at the Hiroshima Prefecture Kure City Dental 

Practitioner's Association made diagnosis using 100 copies of the same panoramic 
radiograph to measure diagnostic accuracy in the same way. 

The results are- approximately 77% sensitivity, approximately 40% specificity, and 
approximately 58% correct diagnosis rate. This allows further accurate screening of a 
10 person at risk of osteoporosis than by practicing dental practitioners. 

Note that the present invention is not limited to the above embodiment, and 
various modifications are possible within the scope of the present invention. 



Industrial Applicability 

15 Since a necessary operation for diagnosis can be reduced to being only specifying a 

body part with the osteoporosis diagnosis aiding apparatus utilizing panoramic 
radiographs according to the present invention, diagnostic ability of the examiner is not 
particularly questioned. 

Furthermore, implementing this apparatus in digital imaging devices being 

20 popularized will allow wide usage thereof at private dental clinics. 

Panoramic radiography is often used for odontopathy diagnosis at dental clinics 
and is practiced with higher frequency than other types of radiography used for medical 
use. Therefore, more persons without concern for osteoporosis may be screened as 
compared to using the conventional methods. 

25 After screening with the apparatus according to the present invention, patients 

with a risk of osteoporosis are instructed to undergo DXA testing at another facility, 
thereby contributing to early detection and early treatment. 
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